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Abstract. We investigate the equilibrium and nonequilibrium properties of the core 
region of vortices in layered superconductors. We discuss the electronic structure of 
singly and doubly quantized vortices for both s-wave and d-wave pairing symmetry. 
We consider the intermediate clean regime, where the vortex-core bound states are 
broadened into resonances with a width comparable to or larger than the quantized 
energy level spacing, and calculate the response of a vortex core to an a.c. electro- 
magnetic field for vortices that are pinned to a metallic defect. We concentrate on the 
case where the vortex motion is nonstationary and can be treated by linear response 
theory. The response of the order parameter, impurity self energy, induced fields and 
currents are obtained by a self-consistent calculation of the distribution functions and 
the excitation spectrum. We then obtain the dynamical conductivity, spatially resolved 
in the region of the core, for external frequencies in the range, O.IZi < hui ~^3A. We 
also calculate the dynamically induced charge distribution in the vicinity of the core. 
This charge density is related to the nonequilibrium response of the bound states and 
collective mode, and dominates the electromagnetic response of the vortex core. 

1 Introduction 

Vortex motion is the principle mechanism for resistive losses in type II supercon- 
ductors. Vortices also provide valuable information about the nature of low lying 
excitations in the superconducting state. In clean s-wave BCS superconductors 
the low-lying excitations in the core are the bound states of Caroli, de Gennes 
and Matricon jlj. These excitations have superconducting as well as normal 
metallic properties. For example, these states are the source of circulating su- 
percurrents in the equilibrium vortex core, and they are strongly coupled to the 
condensate by Andreev scattering Furthermore, the response of the vortex 
core states to an electromagnetic field is generally very different from that of nor- 
mal electrons. However, in the dirty limit, H/t 3> A, the the Bardeen-Stephen 
model Q of a normal-metal spectrum with the local Drude conductivity in the 
core provides a reasonable description of the dissipative dynamics of the vortex 
core. The opposite extreme is the "superclean limit", H/t <§; A 2 /Ef, in which 
the quantization of the vortex-core bound states must be taken into account. In 
this limit a single impurity and its interaction with the vortex core states must 
be considered. The a.c. electromagnetic response is then controlled by selection 
rules governing transition matrix elements for the quantized core levels and the 
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level structure of the core states in the presence of an impurity [HqQ- In the 
case of d-wave superconductors in the superclean limit, nodes in the spectrum 
of bound states lead to a finite dissipation from Landau damping for T — ► || . 

The superclean limit is difficult to achieve even for short coherence length 
superconductors; weak disorder broadens the vortex core levels into a quasi- 
continuum. We investigate the intermediate-clean regime, A 2 /Ef <C H/r <C A, 
where the discrete level structure of the vortex-core states is broadened and the 
selection rules are broken due to strong overlap between the bound state wave 
functions. However, the vortex core states remain well defined on the scale of 
the superconducting gap, A. In this regime we can take advantage of the power 
of the quasiclassical theory of nonequilibrium superconductivity ||9|, |lO| , pl]JT^Jl^] . 

The energy required to maintain a net charge density of order an elementary 
charge per particle within a coherence volume (or coherence area in two dimen- 
sions) is much larger than the condensation energy. Thus, charge accumulation 
in the vortex core is strongly suppressed. In order to reduce the Coulomb energy 
associated with the charge accumulation an internal electrochemical potential, 
<P(R; t), develops in response to an external electric field. This potential produces 
an internal electric field, E iat (R; t), which is the same order of magnitude as the 
external field. Even though the external field may vary on a scale that is large 
compared to the coherence length, £o> the internal field develops on the coherence 
length scale. The source of the internal field is a charge density that accumulates 
inhomogeneously over length scales of order the coherence length. It is necessary 
to calculate the induced potential self consistently from the spatially varying 
order parameter, spectral function and distribution function for the electronic 
states in the vicinity of the vortex core. An order of magnitude estimate shows 
that to produce an induced field of the order of the external field, the dynami- 
cally induced charge is of order e (A/Ef)(5v UJ /A), where Sv^ ~ eE° xt / £ uj is the 
typical energy scale set by the strength of the external field. This charge density 
accumulates predominantly in the vortex core region and creates a dipolar field 
around the vortex core. For a pinned vortex the charge accumulates near the 
interface separating the metallic inclusion from the superconductor. 

Disorder plays a central role in the dissipative dynamics of the mixed state 
of type II superconductors. Impurities and defects are a source of scattering that 
limits the mean free path of carriers, thus increasing the resistivity. Defects also 
provide 'pinning sites' that inhibit vortex motion and suppresses the flux-flow re- 
sistivity. However, for a.c. fields even pinned vortices are sources for dissipation. 
The magnitude and frequency dependence of this dissipation depends on the 
electronic structure and dynamics of the core states of the pinned vortex. In the 
analysis presented below we consider vortices in the presence of pinning centers. 
We model a pinning center as a normal metallic inclusion which is coupled to the 
electronic states of the superconductor through a highly transmitting interface. 
In this model the charge dynamics of the electronic states near the interface 
between the pinning region and superconductor plays an important role in the 
electromagnetic response of the core. 
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In the next section we provide a short summary of the noncquilibrium quasi- 
classical equations, including the transport equations for the quasiparticle distri- 
bution and spectral functions, constitutive equations for the order parameter, im- 
purity self-energy and electromagnetic potentials. In section || we present calcu- 
lations for the the electronic structure of vortices for superconductors with both 
s-wave and d-wave pairing symmetry. The results are based on self-consistent 
calculations of the order parameter and impurity self energy. For s-wave su- 
perconductors, impurity scattering leads to inhomogeneous broadening of the 
vortex core bound states, as well as bands of impurity states within the gap. 
In the case of <i-wave pairing the core states are further broadened by coupling 
between bound states and the continuum states through impurity scattering. 
We also discuss the structure of doubly quantized vortices and vortices bound to 
mesoscopic size metallic inclusions. In the case of the doubly quantized vortex 
there are two branches of zero-energy bound states centered at finite impact 
parameter from the vortex center. This leads to a unique signature of a doubly 
quantized vortex: currents in the core circulate opposite to the supercurrents 
outside the core. Section ^ summarizes calculations of the vortex core dynam- 
ics for s-wave vortices in the presence of impurity scattering. We describe the 
charge dynamics of the vortex core for both pinned and unpinned vortices, and 
calculate the local a.c. conductivity that results from the coupled dynamics of 
the order parameter collective mode and the quasiparticle bound states in the 
vortex core. We discuss energy transport by the core states and the absorption 
features in the conductivity spectrum, which we interpret in terms of absorption 
within the bound-state band centered at the Fermi level and resonant transitions 
involving the bound and continuum states. 

2 Nonequilibrium Transport Equations 

The quasiclassical theory describes equilibrium and noncquilibrium properties of 
superconductors on length scales that are large compared to microscopic scales 
(i.e. the lattice constant, Fermi wavelength, kj 1 , Thomas-Fermi screening length, 
etc.) and energies that are small compared to the atomic scales (e.g. Fermi en- 
ergy, Ef, plasma frequency, conduction band width, etc.). Thus, there are small 
dimensionless parameters that define the limits of validity of the quasiclassical 
theory In particular, we require fc/£o 3> 1, ksTc/Ef <C 1 and hut <C Ef, where 
the a.c. frequencies of interest are typically of order A ~ T c , or smaller, and the 
length scales of interest are of order the coherence length, £o = hvf /2TtkBT c , 
or longer. Hereafter we use units in which % = ks = 1, and adopt the sign 
convention e = — \e\ for the electron charge. 

In quasiclassical theory quasiparticle wavepackets move along nearly straight, 
classical trajectories at the Fermi velocity. The classical dynamics of the quasi- 
particle excitations is governed by semi-classical transport equations for their 
phase-space distribution function. The quantum mechanical degrees of freedom 
are the "spin" and "particle-hole degree of freedom" , described by 4 x 4 density 
matrices (Nambu matrices). The quantum dynamics is coupled to the classical 
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dynamics of the quasiparticles in phase space through the matrix structure of 
the quasiclassical transport equations. 

The nonequilibrium quasiclassical transport equations ]j||lO 1 1 l^,[l3| are for- 



mulated in terms of a quasiclassical Nambu-Keldysh propagator g(pf,R;e,t), 
which is a matrix in the combined Nambu-Keldysh space, and is a function of 
position _R, time t, energy e, and momenta pf on the Fermi surface.^ We denote 
Nambu-Keldysh matrices by a "check", and their 4x4 Nambu submatrices of 
advanced (A), retarded (R) and Keldysh-type (K) propagators by a "hat". The 
Nambu-Keldysh matrices for the quasiclassical propagator and self-energy have 
the form, 

4« *k\ t<jR a K 



where g R > A > K are the retarded (R), advanced (A) and Keldysh (K) quasiclassical 
propagators, and similarly for the self-energy functions. Each of these compo- 
nents of g and a are 4x4 Nambu matrices in combined particle-hole-spin space. 
For a review of the methods and an introduction to the notation we refer to Refs. 
[ fn5||l4| , |l6f . In the compact Nambu-Keldysh notation the transport equations and 
the normalization conditions read 

[(e + -v f ■ A)f 3 - eZ <Pl - A mf - v mf - a, ,g]^ + iv f ■ Vg = , (2) 

g®g = -n 2 l , (3) 
where the commutator is [A, B]$ = A® B — B ® A, 

I®B(M) = e5( 8 ' 9 '- 8 ' 8 '»i(e,i)B(e,t). (4) 

The vector potential, A(R;t), includes A (R) which generates the static mag- 
netic field, Bq(R) = V x Aq(R), as well as the non-stationary vector potential 
describing the time- varying electromagnetic field; A m t(pf,R;t) is the mean- 
field order parameter matrix, P m f(pf,R;t) describes diagonal mean fields due 
to quasiparticle interactions (Landau interactions), and ai(j)f, R; e, t) is the im- 
purity self-energy The electrochemical potential ^(.R; t) includes the field gen- 
erated by the induced charge density, p(R;t). The coupling of quasiparticles 
to the external potential involves virtual high-energy processes, which result 
from polarization of the non-quasiparticle background. The interaction of quasi- 
particles with both the external potential 4> and the polarized background can 
be described by coupling to an effective potential Zq<P Q). The high-energy 
renormalization factor Zq is defined below in Eq. ([12| ). The coupling of the 
quasiparticle current to the vector potential in Eq. (0) is given in terms of the 
quasiparticle Fermi velocity. No additional renormalization is needed to account 
for the effective coupling of the charge current to the vector potential because 
the renormalization by the non-quasiparticle background is accounted for by 

1 In quasiclassical theory the description in terms of the variables (e,py,R) is related 
to the phase-space description in p — R space by a transformation, <?(p/ , e; R, t) = 
/(p, R; t), with e — e(p, R; t) — /x and p = pf [jpjj]. 
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the effective potentials that determine the band structure, and therefore the 
quasiparticle Fermi velocity. 

2.1 Constitutive Equations 

Equations must be supplemented by Maxwell's equations for the electro- 

magnetic potentials, and by self-consistency equations for the order parameter 
and the impurity self-energy. We use the weak-coupling gap equation to describe 
the superconducting state, including unconventional pairing. The mean field self 
energies are then given by, 

A%f(p f ,R;t) = N f [ +eC ^-(V(pf,p f )f K (p f ,R;e,t)) , (5) 

0*f(p f ,R;t) = N f ^-(A(p f ,pi f )g K (p> f ,R;e,t)) , (6) 

AK f (pf,R;t) = , v% f (p f ,R;t) = . (7) 
The impurity self-energy, 

&i(pf,R;e,t) = n t t(pf,p f ,R;e,t) , (8) 

is specified by the impurity concentration, rii, and impurity scattering t- matrix, 
which is obtained from the the self-consistent solution of the i-matrix equations, 

i(Pf,p'},R;e,t) = u(jpf,p" f ) 

+ Nf(u(pf,Pf)®g( P f,R;e,t)®t( Pf ,p'},R;e,t)) . (9) 

The Nambu matrix f K is the off-diagonal part of g K , while g K is the diagonal 
part in particle-hole space. The Fermi surface average is defined by 

/ If d 2 P'f f d 2 P't 

<•■•>- wJwfwi ( -" ) ' N,= l^m-\' m 

where Nf is the average density of states on the Fermi surface. The other 
material parameters that enter the self-consistency equations are the dimen- 
sionless pairing interaction, N fV(pf ,p'f), the dimensionless Landau interaction, 
NfA(pf,p'j), the impurity concentration, rij, the impurity potential, u(pf,p'f), 
and the Fermi surface data: pf (Fermi surface), Vf(pf) (Fermi velocity). We 
eliminate both the magnitude of the pairing interaction and the cut-off, e c , in 
favor of the transition temperature, T c , using the linearized, equilibrium form of 
the mean- field gap equation (Eq. (0)). 

The quasiclassical equations are supplemented by constitutive equations for 
the charge density, the current density and the induced electromagnetic poten- 
tials. The formal result for the non-equilibrium charge density, to linear order in 
A/Ef, is given in terms of the Keldysh propagator by 

pW(R;t) = eNf ^-.(Z(p' f ) Tr [g K (p' f ,R;e,t)] > -2e 2 N f Z ^(R;t) , (11) 
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with the renormalization factors given by 

Z(p f ) = l-(A(p' f , Pf )), Z = (Z(p f )). (12) 

The high-energy renormalization factor is related to an average of the scattering 
amplitude on the Fermi surface by a Ward identity that follows from the con- 
servation law for charge filf] . The charge current induced by A(R; t), calculated 
to leading order in A/Ef, is also obtained from the Keldysh propagator, 

j^(R;t) = eNf J ^-Jr(v f (p' f )f 3 g K (p' f ,R;e,t)) . (13) 

There is no additional high-energy renormalization of the coupling to the vector 
potential because the quasiparticle Fermi velocity already includes the high- 
energy renormalization of the charge-current coupling in Eq. [l3| Furthermore, 
the self-consistent solution of the quasiclassical equations for g K ensures the 
continuity equation for charge conservation, 

d tP W(R;t) + V-jM(R;t)=0, (14) 

is satisfied. 

An estimate of the contribution to the charge density from the integral in 
Eq. ( pT| ) leads to the condition of "local charge neutrality" |l^,|l8). A charge 
density given by the elementary charge times the number of states within an 
energy interval A around the Fermi surface implies p^ ~ 2eNfA. Such a charge 
density cannot be maintained within a coherence volume because of the cost in 
Coulomb energy. The Coulomb energy is suppressed by requiring the leading 
order contribution to the charge density vanish: i.e. p^(R;t) = 0. Thus, the 
spatially varying renormalizcd electro-chemical potential, Zq<P, is determined by 

2eZ Q <P(R; t) = J & ° ^-Tr(Z(p' f )g K (p> f , R; e, t)) . (15) 
The continuity equation implies V ■ j^(R; t) = 0. We discuss violations of the 



charge neutrality condition ([15|), which are higher order in A/Ef, in Sec. 4.1 
Finally, Ampere's equation, with the current given by Eq. (|l3|), determines the 
vector potential in the quasiclassical approximation, 

VxVxA(fl;i) = J ^-Tr(v f (p' f )f 3 g K (p' f ,R;e,t)) . (16) 

Equations (|2|)-([)|) and (|l5|)-(|l6|) constitute a complete set of equations for 
calculating the electromagnetic response of vortices in the quasiclassical limit. 
For high-K superconductors we can simplify the self-consistency calculations to 
some degree. Since quasiparticles couple to the vector potential via |«/ - A, Eq. 
( p^| ) shows that this quantity is of order 87re 2 N fv'j / c 2 = 1/A 2 , where A is the 
magnetic penetration depth. Thus, for k = A/£o 3> 1, as in the layered cuprates, 
the feedback effect of the current density on the vector potential is small by 
factor 1/k 2 . 
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2.2 Linear Response 

For sufficiently weak fields we can calculate the electromagnetic response to 
linear order in the external field. The propagator and the self-energies are sep- 
arated into unperturbed equilibrium parts and terms that are first-order in the 
perturbation, 

g = go + 8g , A mf = A + SA mf , a { = <x + Sat , (17) 

and similarly for the electromagnetic potentials, A = Aq + 5A, <P = (5<£. The 
equilibrium propagators obey the matrix transport equation, 

[(e + -v f ■ A )f 3 - A -a ,go]+iVf ■ Vg = . (18) 

These equations are supplemented by the self-consistency equations for the mean 
fields, Eqs. the impurity self energy, Eqs. (^|)-(p|), the local charge- 

neutrality condition for the scalar potential, Eq. (fig), Ampere's equation for 
the vector potential, Eq. ([l6|), the equilibrium normalization conditions, 

.9o = ~^ 2 i - (19) 

and the equilibrium relation between the Keldysh function and equilibrium spec- 
tral density, 

5o*=tanh(^) [g$-g£]. (20) 

The first-order correction to the matrix propagator obeys the linearized trans- 
port equation, 



[(e+-VfAo)f3-Ao-a ,5g] +iVfVSg= [SA mf +5ai+6i> , g }® , (21) 



with source terms on the right-hand side from both the external field (Sv) and 
the internal fields (SA m f, Sdi). In addition, the first-order propagator satisfies 
the "orthogonality condition" , 

5o ® Sg + Sg <g> g = . (22) 

obtained from linearizing the full normalization conditional The system of linear 
equations are supplemented by the equilibrium and first-order self-consistency 
conditions for the order parameter, 

A«' A (p f ,R) = N f ^-(V(p f ,pi f )f K (p' f ,R;e)) , (23) 

6A%f(p f ,R;t) = N f ^-(V(p f ,pi f )Sf K (pi f ,R;e,t)), (24) 



2 Note that the convolution product between an equilibrium and a nonequilibrium 
quantity simplifies after Fourier transforming t — - > u: Ao CS> 5B(e,ui) = Ao(e + 
w/2)B(e,w), B(e,w) <8> Aq = B(e, iv)A (e - w/2). 
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and the impurity self-energy, 

oo (Pf , R, e) = rii t (p f ,p f ,R;e), (25) 
t (Pf,Pf, R, e) = u(pf,p'}) + N f (u(p f ,p' f )g (p' f ,R; e)t (p' f ,Pf, R, e)> , (26) 
Sa z (p f ,R; e, t) = n z N f (i Q (p f ,p' f ,R; e) ® Sg(p' f ,R; e, t) $ t (p' f ,Pf, R, e)) -(27) 

In general the diagonal mean fields also contribute to the response. However, 
we do not expect Landau interactions to lead to qualitatively new phenomena 
for the vortex dynamics, so we have neglected these interactions in the following 
analysis and set A(pf,p'j) — (i.e. v m f = 0). As a result the local charge 
neutrality condition for the electro-chemical potential becomes, 

2e6<P(R; t) = / " ^Tr (6 g K (pf,R; e, t)) . (28) 

In what follows we work in a gauge in which the induced electric field, 
E iad (R;t), is obtained from S ( l ) (R;t) and the uniform external electric field, 
E e ^(t), is determined by the vector potential SA u (t). For A/£o ^ 1 we can 
safely neglect corrections to the vector potential due to the induced current. 
Thus, in the Nambu-Keldysh matrix notation the electromagnetic coupling to 
the quasiparticles is given by 

5v = --v r 5A u (t)f 3 + e5$(R;t)l . (29) 

The validity of linear response theory requires the external perturbation 5v be 
sufficiently small and that the induced vortex motion responds to the external 
field at the frequency set by the external field. At very low frequencies frictional 
damping of the vortex motion, arising from the finite mean free path of quasipar- 
ticles scattering from impurities, gives rise to a nonlinear regime in the dynamical 
response of a vortex. This regime is discussed extensively in the literature 
and is not subject of our study. However, for sufficiently small field strengths the 
vortex motion is nonstationary over any time interval, although it may be re- 
garded as quasi-stationary at low enough frequencies. The nonstationary motion 
of the vortex can be described by linear response theory if Sv <C 1/t for lo 
and 8v -C to for lo Note that the frequency of the perturbation, u, is not 

required to be small compared to the gap frequency; it is only restricted to be 
small compared to atomic scale frequencies, e.g to <C Ef/h. 



Self-consistent solutions of Eqs. (p4|), ( |27[ ) and (|28| ) for the self-energies and 
scalar potential are fundamental to obtaining a physically sensible solution for 
the electromagnetic response. The dynamical self-energy corrections are equiv- 
alent to "vertex corrections" in the Kubo formulation of linear response theory. 
They are particularly important in the context of noncquilibrium phenomena in 
inhomogeneous superconductors.^ In our case these corrections are of vital im- 
portance; the self-consistency conditions enforce charge conservation. In partic- 
ular, Eqs. (p5|)-(p7j) imply charge conservation in scattering processes, whereas 

3 Vertex corrections usually vanish in homogeneous superconductors because of trans- 
lational and rotational symmetries. Inhomogeneous states break these symmetries 
and typically generate non-vanishing vertex corrections. 
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( p3[ ) and ( p4| ) imply charge conservation in particle-hole conversion processes; 
any charge which is lost (gained) in a particle-hole conversion process is com- 
pensated by a corresponding gain (loss) of condensate charge. It is the coupled 
quasiparticle and condensate dynamics which conserves charge in superconduc- 
tors. Neglecting the dynamics of cither component, or using a non-conserving 
approximation for the coupling leads to unphysical results. 

Self-consistent calculations for the equilibrium order parameter, impurity 
self-energy and local excitation spectrum (spectral density) are necessary inputs 
to the linearized transport equations for the dynamical response of a vortex. The 
equilibrium spectral function also provides key information for the interpreta- 
tion of the dynamical response. Because of particle-hole coherence the spectral 
density is sensitive to the phase winding and symmetry of the order parameter, 
as well as material properties such as the transport mean-free path and impu- 
rity cross-section. In the following section we present results for the low-energy 
excitation spectra of singly- and doubly-quantized vortices in layered supercon- 
ductors with s-wave and d-wave pairing symmetry. 



3 Electronic Structure of Vortices 

The local density of states for excitations with Fermi momentum pt is obtained 
from the retarded and advanced quasiclassical propagators, 

N(p f ,R; e) = N f ^-Tr [hg£(p f , R e) - T 3 g«(p f: R; e)] . (30) 

This function measures the local density of the quasiparticle states with energy 
e at the point P f on the Fermi surface. The local density of states (LDOS) is 
obtained by averaging this quantity over all momentum directions of the quasi- 
particles, 

N(R;e) = (N(p' f ,R;e)) . (31) 

The product of the angle-resolved density of states and the Fermi velocity, Vf, 
determines the current density carried by these states; Vf also defines the di- 
rection of a quasiclassical trajectory passing through the space point R. We 
introduce the angle-resolved spectral current density || , 

j(p f ,R;e) = 2ev f (p f )N{p f ,TL;e), (32) 

which measures the current density carried by quasiparticle states with energy 
e moving along the trajectory defined by Vf. The local spectral current density 
is then defined as 

j(R;e) = (j(p' f ,R;e)) , (33) 

and the total current density is obtained by summing over the occupied states 
for each trajectory, 

i(R) = J def(e)j(R;e), (34) 
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where /(e) = 1/(1 + e' 36 ). Self-consistent calculations of the equilibrium structure 
and spectral properties of vortices are relatively straight- forward computations. 
Below we present results for s-wave and d-wavc pairing symmetry with impurity 
scattering included. 

The calculations reported are carried for a circular Fermi surface, with an 
isotropic Fermi momentum pf and Fermi velocity Vf. The elastic scattering rate 
is chosen to represent the intermediate-clean regime, A 2 /Ef < h/r <C A. The 
pairing potential can be represented as a sum over invariant products of basis 
functions \r\i \i(p : f)\i — 1 . . . dr} for the irreducible representations of the crystal 
point group labeled by T, 

N f V( Pf ,p' f ) = VrVh(Pf)vMp' f ) ■ (35) 

The pairing interaction, vr, and the cutoff, e c , are eliminated in favor of the 
instability temperature for pairing in symmetry channel r. We limit the discus- 
sion here to even-parity, one-dimensional representations, which for tetragonal 
symmetry includes the 's-wave' (identity) representation, Ai s , two 'd-wave' rep- 
resentations, Big and B2 g , and a 'g-wave' representation, A2 g . The corresponding 
basis functions we use are listed in Table 

Table 1. Symmetry classes and model basis functions for the ID 
even-parity representations of D^h- The angle <f) is the angular 
position of p/ on the Fermi surface with respect to the crystallo- 
graphic a-axis (=a>axis). 



Pairing 


Symmetry 


Representation [F] 


Basis Function [r/r] 


s- 


wave 


Aig 


1 


d- 


wave 


Big 


V / 2cos(20) 


d' 


-wave 


B 2g 


v^sin^) 


9- 


wave 


A 2g 


\/2sm(4:(j>) 



The results for the order parameter, impurity self energy and spectral prop- 
erties of vortices that follow are calculated self consistently in the i-matrix ap- 
proximation for point impurities (pure s-wave scattering), i.e. u(pf : p'j) = u^l. 
The quasiparticle scattering rate, l/2r, and normalized impurity cross section, 
a, are then given by, 

— -^-n _ (irN f u ) 2 

2r nN f ' 1 + (irN f uo) 2 ' { ' 

The Born limit corresponds to a <C 1, while unitary scattering corresponds to 
| i*o | — > oo or a — > 1. In the calculations that follow the temperature is set at 
T = 0.3T C , and the mean free path is is chosen to represent the intermediate- 
clean regime; I = 10 £o- 
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3.1 Singly quantized vortices for S-wave Pairing 

For isotropic s-wave pairing the equilibrium order parameter for an isolated 
vortex with winding number p has the form, 

A(p f ,R) = \A(R)\e i ™, (37) 

where the amplitude |Z\(i?)| is isotropic and if is the azimuthal angle of R in 
the plane. 




-A A 8 -A A e 

Fig. 1. Angle-resolved density of states for an s-wave vortex for quasiparticles propa- 
gating along trajectories parallel to the :r-direction, at points along the j/-axis spaced 
by 1.38 £o- The left panel is for Born scattering, and the right panel is for unitary 
scattering with the same mean free path of £ — 10£o- The temperature is T = 0.3T C . 

The angle-resolved local density of states spectra for a singly quantized vortex 
is shown in Fig. [I] for space points R = (0, y) along the y-axis with a spacing of 
Sy = \/3n/4i;o ~ 1.36£o, and for trajectories parallel to the x-axis, Vf = vje x . 
The vortex center is at y — 0, and the phase winding is such that the direction 
of the superfiow is in x-direction for spectra with negative y coordinate. 

Far outside the core the angle-resolved density of states resembles the BCS 
density of states with a gap in the spectrum roughly between e = ±Z\, and peaks 
in the spectrum near the continuum edges. Careful inspection of Fig. [I] shows 
that the coherence peak for positive energy at y = —&5y ~ — H£o 1S n °t at 
e = A, but is shifted to higher energy by the Doppler effect, Ac = v/ • p s , where 
p s = |fi,V$— - A is the condensate momentum (3). In a homogeneous superfiow 
field the spectrum is the Doppler-shiftcd BCS spectrum; the Doppler shift in- 
creases with condensate momentum until pairbreaking sets in at the bulk critical 
momentum, p c — A/vt. However, in the vortex core nonlocal effects associated 
with the inhomogeneous flow field lead to a redistribution of the spectral weight 
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near the gap edge. The positive energy continuum edge is broadened consid- 
erably compared to the square-root singularity in the absence of the Doppler 
effect. The continuum starts at +A even for the Doppler-shifted spectra near 
the maximal current regions. In contrast, the negative energy continuum edge 
shows sharp structures due to the accumulation of spectral weight in the region 
between —A and — A + vy • p s . The sharp structure corresponds to a bound 
state that is separated from the continuum edge. The density of states at the 
continuum edge drops precipitously at e = — A. We emphasize that nonlocal ef- 
fects lead to qualitative differences in the spectrum near the gap edges compared 
to the widely used approximation of Doppler-shifted quasiparticles in a locally 
homogeneous superfiow fl9|| . In clean superconductors nonlocal effects dominate 
the spectrum. 

For a homogeneous superfiow the current is carried mainly by the states 
that are Doppler shifted the region between —A and A. The spectral current 
density shows that contributions to the current density from states outside of 
this region nearly cancel. In the case of a vortex, the bound state that splits 
off from the continuum not only robs the Doppler-shifted continuum edge of its 
spectral weight, but the bound state also carries most of the supercurrent in 
the vortex core region. At distances approaching the vortex center the bound 
state is clearly resolved and disperses through e = at zero impact parameter. 
As shown in the left panel of Fig. [I] the bound state also broadens considerably 
for Born scattering as it disperses towards the Fermi level, but remains a sharp 
resonance in the limit of unitary scattering, as shown in Fig. Ill 

The coherence peaks are completely suppressed at the vortex center, both in 
the Born and unitary limits. Note the difference in the evolution of the spectral 
weight for the bound and continuum states: for positive energies the spectral 
weight of the coherence peak is shifted to higher energies as one approaches the 
vortex core, with a continuously diminishing intensity at the coherence peak. In 
contrast, the spectral weight of the negative energy coherence peak is transferred 
to the bound state which splits off from the continuum. In principle there could 
be additional, secondary bound states, which would split off from the continuum 
if the vortex core were wider. However, a self-consistent calculation of the order 
parameter suppresses the secondary bound state in the vortex core. We return 
to the discussion of secondary bound states when we discuss the spectrum of 
pinned vortices. For trajectories with y > the structure of the spectrum for 
positive and negative energies is reversed because the superfiow is now counter 
to the quasiparticle velocity, leading to negative Doppler shifts and bound states 
at positive energies below the continuum edge. The small spectral features that 
appear at the energies corresponding to the negative of the bound state energies 
are due to mixing of with trajectories of opposite velocity by backscattering from 
impurities. 

Finally, consider the differences in the spectral features for unitary versus 
Born scattering. In Fig. || we show the local density of states at the center of the 
vortex for several scattering cross sections ranging from the Born limit (ct < 1) 
to the unitary limit [d = 1). In the Born limit the broadening is a maximum. 
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Fig. 2. Local density of states (LDOS) in the center of an s-wave vortex for different 
effective impurity scattering cross sections a from Born to unitary limit, with a fixed 
mean free path I — 10£o, and temperature T — 0.3T C . From top to bottom: a = 
0.0,0.2,0.4,0.6,0.8,1.0. 



For small, but finite cross section two bands of impurity bound states split off 
from the zero-energy resonance and remove spectral weight from the central 
peak. The impurity bands evolve towards the continuum edges as the cross- 
section increases and merge with them in the unitary limit. When the impurity 
bound-state bands no longer overlap the central peak the zero-energy resonance 
sharpens dramatically with the width remaining constant as the unitary limit is 
approached. The overlap between the impurity bands and the continuum edges 
in the unitary limit is determined by the scattering rate, and is increasing with 
increasing scattering rate, 1/t. As shown in the right panel of Fig. El the impurity 
bands are localized in the vortex core region; their existence depends on impurity 
scattering in a region where the phase changes rapidly over length scales of order 
the coherence length. 

3.2 Singly quantized vortices for D-wave Pairing 

For ci-wave pairing symmetry the order parameter has the form, 

A(p f ,R) = r, Blg (p f )A(K), (38) 

where ?y Bl (p/) changes sign on the Fermi surface at the points, Pf x — ip/y 
These nodal points lead to strong anisotropy and gapless excitations in the 
quasiparticle spectrum, which feeds back to produce anisotropy in the spatial 
structure of the order parameter in the core region of of a vortex. The spatial 
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part of the order parameter, 

A(R) = \A(H)\e M(R \ (39) 

for a vortex, at distances far from the core, approaches that of an isotropic vortex: 
|Z\(R)| — > A(T) and i9 — > pip. However, in the core region the current density is 
comparable to the critical current density and develops a four-fold anisotropy as 
a result of the backflow current concentrated near the nodes j2(j] . This current- 
induced pairbreaking effect is dominant for flow parallel to the nodal directions 
and leads to weak anisotropy of the order parameter in the core region. 

The electronic structure of the d-wave vortex, presented in Fig. ^, shows 
distinct differences from that of an s-wave vortex. The resonances in the vortex 
core are broader than the core states of vortex in an s-wave superconductor. 
Mixing with extended states in nodal direction broadens the peaks near the 
continuum edges. Also note the effect of impurity scattering on the spectra far 
from the vortex core region, which show a broadened continuum rather than a 
sharp continuum edge. In the right panel of Fig. [| we show the angle-resolved 
density of states for a <i-wave vortex for unitary impurity scattering. As in the 
s-wave case there is sharpening of the zero resonance, but note that the width 
of the resonance is much broader than in the case of s-wave pairing because 
impurity scattering provides a coupling and mixing of the bound state resonance 
with the low-energy extended states for momenta near the nodal directions, even 
for a trajectory in the antinodal direction. 




Fig. 3. Angle resolved density of states for a d-wave vortex for a quasiparticle trajec- 
tory along the a;-direction (antinodal) as a function of impact parameter (y-direction) . 
The impact parameter spacing is Sy = 1.36£q- The left panel corresponds to impu- 
rity scattering in the Born limit, and the right panel is for the unitary limit. The 
temperature is T = 0.3T C . 



In Fig. |] we show the local density of states (LDOS) for both s-wave (top 
two panels) and d-wave pairing (bottom two panels) as a function of position R 
along an anti-nodal direction (bottom left) and along a nodal direction (bottom 
right). The LDOS is obtained by averaging the angle-resolved density of states 
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over the Fermi surface at a particular space point R. The averaging, together 
with the dispersion of the bound states and resonances as a function of angle 
leads to one-dimensional bands with characteristic Van Hove singularities, which 
are clearly visible for s-wave symmetry, but considerably smeared and broadened 
for d-wave symmetry. The LDOS for a vortex with s-wave pairing symmetry is 
shown for impurity scattering in both the Born (top left) and unitary limits 
(top right). The spectra show the characteristic bound state bands, Van Hove 
singularities and the dramatic reduction in the width of the zero bias resonance 
for unitary scattering. 
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Fig. 4. Local density of states for an s-wave vortex (top panels) and a d-wave vortex 
(bottom panels), as a function of distance from the vortex center with a spacing 5y = 
0.79£o- For s-wave pairing the Born limit is shown on the left and unitary scattering 
is on the right. For d-wave pairing the LDOS is shown only for Born scattering, but 
along axes parallel to a anti-nodal direction (left) and parallel to the nodal (right). The 
temperature is T = 0.3T C . 

Calculations of the LDOS in the superclean limit for a vortex with d-wave 
pairing are discussed by Schopohl and Maki and by Ichioka et al. [p2[ . The 
self-consistent calculations shown in Fig. H include impurity scattering in the 
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Born limit for space points, R, along two different directions; the left panel 
corresponds the LDOS measured as a function of distance along the anti-nodal 
direction and the right panel is the LDOS measured along the nodal direction. 
The nodes of the order parameter for ci-wave pairing lead to continuum states 
with energies down to the Fermi level. These states are visible in the LDOS as 
the smooth background extending to zero energy from both positive and negative 
energies, even for distances far from the core. In the vortex core region several 
broad peaks disperse as a function of distance from the vortex center. These 
peaks correspond to broadened Van Hove singularities resulting from averaging 
the vortex core resonances over the Fermi surface for at a fixed position R. The 
differences in the <i-wave spectra for the two directions reflects the weak fourfold 
anisotropy of the LDOS around the vortex at fixed energy. 

3.3 Vortices pinned to mesoscopic metallic inclusions 

The calculations discussed above describe the average effects of atomic scale im- 
purity disorder on the spectral properties of a vortex. A specific defect can also 
act as a pinning site for a vortex. We model such a defect as a mesoscopic size, 
normal metallic inclusion in the superconductor. The defect is assumed to be a 
circular inclusion, with a radius, £ pin , of order the coherence length of the super- 
conductor. For simplicity we assume that the metallic properties of the inclusion, 
e.g. Fermi surface parameters, are the same as those of the normal-state of the 
host superconductor. The calculations presented below neglect normal reflection 
processes at the interface between the inclusion and the host metal, but include 
Andreev reflection. The analysis and calculations can be generalized for more 
detailed models of a pinning center. 

Figures || and || show the angle-resolved local density of states and the LDOS 
for a vortex pinned on a metallic inclusion of radius £ pin = 7r£o . The order param- 
eter, impurity self energy and spectral densities were calculated self consistently 
for impurity scattering in the Born limit. Qualitative changes resulting from 
the inclusion occur inside the pinning center. The shape of the bound state 
resonance lines are asymmetric in energy. The asymmetry arises from multiple 
Andreev reflection processes by the interface between the pinning center and the 
superconductor, which leads to additional bands of resonances that overlap the 
vortex core resonances. 

In addition to the asymmetry in the linewidth of the resonances the zero- 
energy bound state at the vortex center has a peculiar spectral shape, shown in 
more detail in Fig. [7] for s-wave pairing symmetry, but also visible in the right 
panel of Fig. ^ for d-wave pairing symmetry as well. In contrast to the spectra 
for a vortex without a pinning center, the coherence peaks at the continuum 
edges are present at the vortex center. 

Figure |?j shows the peaks near the continuum edges in the vortex center for 
a pinned vortex. The spectral weight near the continuum edge is taken from 
the resonance at the Fermi level. The enhanced weight at the continuum edges 
is a precursor to the formation of a secondary bound state that splits off from 
the continuum. This can be seen in the evolution of the spectral density at 
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Fig. 5. Angle resolved density of states for a quasiparticle trajectory along the x- 
direction, as a function of impact parameter (y-axis) with a spacing of Sy = 1.36£o 
for a vortex centered on a pinning center of radius of £ pin = 7r£o- The left panel is 
for s-wave symmetry and right panel is for d-wave symmetry. Impurity scattering is 
included for Born scattering with I — 10£o- The temperature is T = 0.3T C . 
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Fig. 6. Local density of states as a function of distance from the vortex center, with a 
spacing of Sy — 0.79£o, for a vortex pinned to an inclusion of radius of £ p = 7r£o- The 
left panel corresponds to s-wave symmetry and the right panel is for d-wave symmetry 
along the anti-nodal direction. The temperature is T = 0.3T C . 
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the vortex center as a function of the radius of the pinning center. There is a 
zero energy resonance for all pinning radii, however, increasing the radius of the 
pinning center transfers spectral weight from the zero-energy resonance to the 
continuum edge. A coherence peak develops, splits off from the continuum edge, 
strengthens and evolves to energies within the gap as the pinning radius changes 
from £ pill = 0.79£o to £ pin = 4.71£o- Thus, the appearance of the coherence peak 
for pinning centers the size of a coherence length or so is a precursor to the 
formation of a secondary bound state within the gap. The spectral weight comes 
at the expense of states just above the continuum edge and the zero-energy 
bound state, which is diminished in intensity with increasing pinning radius. 




Fig. 7. Local density of states at the center of a vortex pinned in the center of a 
normal inclusion. The spectra (from top to bottom) correspond to different pinning 
radii: £ pin = 0, 0.79£o, 1.57&, 2.36£o, 3.14£ , 3.93£ , and 4.71£ - The temperature is 
T = 0.3T C . 



3.4 Doubly quantized vortices 

Vortices with winding numbers larger than p = 1 generally have line energies 
per unit winding number that are larger than that of a singly quantized vor- 
tex.[|] Nevertheless, doubly quantized vortices once formed are metastable and 
in principle it should be possible to observe the rare doubly quantized vortex 
using an atomic probe such as a scanning tunneling microscopic |p5| , p6| . The 
spectrum of a doubly quantized vortex differs in a fundamental way from that 

4 There are counter examples for unconventional pairing with a multi-component order 
parameter in which the lowest energy vortex states are doubly quantized vortices 
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of a singly quantized vortex. The singly quantized vortex, has a single branch of 
states that disperse through zero energy at the vortex center. The zero mode is 
guaranteed in the quasiclassical limit by ir change along trajectories that pass 
through the vortex center. In contrast, there is no phase change along a tra- 
jectory through the vortex center for a doubly quantized vortex, and thus no 
topological requirement enforcing a zero energy bound state at the center of a 
doubly quantized vortex. Nevertheless, there is a spectrum of bound states in 
the cores of doubly quantized vortices which lead to characteristic structures in 
the LDOS and current spectral density of a doubly quantized vortex. 




Fig. 8. Dispersion of the bound states (below the continuum edge) and of the coherence 
peaks (above the continuum edge) for quasiparticle trajectories in x direction as a 
function of y, for a singly quantized vortex (left) and for a doubly quantized vortex 
(right). The center of the vortex is at y = 0. Born impurity scattering with a mean 
free path of i — 10£o is assumed. The temperature is T = 0.3T C . The thin lines show 
the energies of the continuum edges and the Doppler shifted energies of the coherence 
peaks, assuming the London form for the condensate momentum p s = pe^/2r, where 
p is the winding number. The dispersion is shown for s-wave pairing; the corresponding 
data for d-wave pairing is similar. 



This structure was discussed for a doubly quantized vortex in the superclean 
limit in Rcf . || . Figure ^| shows the angle-resolved density of states for a doubly 
quantized vortex in an s-wave superconductor for trajectories parallel to the 
a;-axis at different impact parameters along the y-direction. Two branches of 
vortex bound states cross the Fermi level at distances of order a coherence length 
from the vortex center. Thus, zero-energy bound states exist in the core but 
they are localized (for s-wave pairing symmetry) on trajectories at finite impact 
parameter from the vortex center. The locus of these trajectories forms a circle 
of radius r ba ~ 2.5£o around the vortex center. Also note that the Doppler shift 
of the continuum spectrum is twice that for singly quantized vortices. 

This can be seen by comparing the spectra far from the core in Figs. ^ and 
[l]. The distinctive features in the spectrum of a doubly quantized vortex are 
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Fig. 9. Angle resolved density of states for a doubly quantized vortex along quasipar- 
ticle trajectories parallel to the ^-direction, as a function of impact parameter, y, with 
a spacing Sy = 1.36£o- The left panel corresponds to s-wave symmetry and the right 
panel is for d-wave symmetry. Impurity scattering is included self-consistently in the 
Born limit with a mean free path £ — 10£o- The temperature is T — 0.3T C . 



present for either s- and <i-wave pairing symmetry. The doubling of bound state 
branches that cross the Fermi level is observed for both pairing symmetries. The 
main difference in the spectra, as in the case of singly quantized vortices, is in 
the width of the resonances. 

The most distinguishing feature of doubly quantized vortices is that the su- 
percurrents near the center of the vortex flow counter to the asymptotic superflow 
associated with the phase winding around the vortex [0. The counter circulating 
currents in the core, shown in the right panel of Fig. are due to the bound 
states interior to the radius defined by the zero-energy bound state (Note that 
the zero-energy bound state itself does not carry current). This structure is re- 
vealed in Figs. |l^ and [Oj which show the LDOS and spectral current density 
for a doubly quantized vortex with s-wave and d-wave pairing symmetry. The 
left panel of Fig. |ll] also shows the cumulative spectral current density and the 
reversal of the current for as one branch of bound states disperses below the 
Fermi level for r < r bs . 

These spectra show that the current density near the vortex core is carried 
mainly by the bound states, and that the reversal of the current direction near 
the vortex center is due to the branch of counter-flowing bound states dispersing 
below the Fermi level for impact parameter r < r b5 . In this region of the core the 
co-moving bound state is above the Fermi level so the asymmetry in the level 
occupation produces a counter flowing current. 
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Fig. 10. Local density of states for doubly quantized vortices along trajectories parallel 
to the x-axis as a function of impact parameter along the y-axis with a spacing of 
5y = 1.36£o- The left panel corresponds to s-wave pairing symmetry, and the right panel 
is for d-wave pairing along the anti-nodal direction. Impurity scattering is included in 
the Born limit with a mean free path of £ = 10£o and the temperature is T = 0.3T C . 
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Fig. 11. The spectral current density for a doubly quantized vortex is shown along 
a trajectory parallel to the a;-axis as a function of impact parameter (j/-axis) with a 
spacing Sy — 1.36£o, starting at y = 1.36£o- The thick curve is the cumulative spectral 
current density as a function of e, obtained by integrating y) from — oo to e. The 
right panel is the x-component of the current density as a function of impact parameter 
for s-wave (solid line) and d-wave (dashed line) pairing. 
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4 Nonequilibrium Response 

The dynamics of the electronic excitations of the vortex core play a key role in 
the dissipative processes in type II superconductors. Except in the dirty limit, 
I <C £o , the response of the core states to an electromagnetic field is generally 
very different from that of normal electrons. It is energetically unfavorable to 
maintain a charge density of the order of an elementary charge over a region with 
diameter of order the coherence length. Instead an electrochemical potential 
is induced which ensures that almost no net charge accumulates in the core 
region. However, a dipolar-like charge distribution develops which generates an 
internal electric field in the core. The internal field varies on the scale of the 
coherence length. This leads to a nonlocal response of the quasiparticles to the 
total electric field, even when the applied field varies on a much longer length 
scale and can considered to be homogeneous. The dynamical response of the 
vortex core includes the collective mode of the inhomogeneous order parameter. 
This mode couples to the electro-chemical potential, S^, in the vortex core region. 
This potential is generated by the charge dynamics of vortex core states and 
gives rise to internal electric fields which in turn drive the current density and 
the order parameter near the vortex core region. The induced electric fields in 
the core are the same order of magnitude as the external field. The dynamics of 
the core states are strongly coupled to the charge current and collective mode of 
the order parameter. Thus, the determination of the induced order parameter, 
as well as the spectrum and distribution function for the core states and non- 
equilibrium impurity scattering processes requires dynamical self consistency. 
Numerous calculations of the a.c. response neglect the self-consistent coupling 
of the collective mode and the spectral dynamics, or concentrate on the u) — ► 
3(J . Quasiclassical theory is the only formulation of the theory of 
nonequilibrium superconductivity presently capable of describing the nonlocal 
response of the order parameter and quasiparticle dynamics in the presence of 
mesoscopic inhomogeneities and disorder. The numerical solution to the self- 
consistency problem was presented for unpinned vortices in Ref. [ BlJ . Here we 
report results for the the electromagnetic response of isolated vortices bound to 
a pinning center in a superconductor with s-wave pairing symmetry. 

4.1 Dynamical charge response 

The charge dynamics of layered superconductors has two distinct origins. The 
c-axis dynamics is determined by the Josephson coupling between the conduct- 
ing planes. Here we are concerned with the in-plane electrodynamics associated 
with the response of the order parameter and quasiparticle states bound to the 
vortex core. We assume strong Josephson coupling between different layers, and 
neglect variations of the response between different layers. This requires that 
the polarization of the electric field be in-plane, so that there is no coupling of 
the in-plane dynamics to the Josephson plasma modes. The external electromag- 
netic field is assumed to be long wavelength compared to the size of the vortex 
core, A EM 3> £o- In this limit we can assume the a.c. electric field to be uniform 
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and described by a vector potential, E w (t) = — ~dtA w . We can also neglect the 
response to the a.c. magnetic field in the limit A>(o, In this case the spatial 
variation of the induced electric field occurs mainly within each conducting layer 
on the scale of the coherence length, £o- Poisson's equation implies that induced 
charge densities are of order S<P/£,q, where 5<1> is the induced electrochemical po- 
tential in the core. This leads to a dynamical charge of order e (A/Ef) in the 
vortex core. Once the electrochemical potential is calculated from Eq. (^8|) we 
can calculate the charge density fluctuations of order (A/Ef) 3 from Poisson's 
equation, 

P ^(R;t) = -±V 2 <P(R;t). (40) 




Fig. 12. Left: Total electric field in the center of a pinned s-wave vortex as a func- 
tion of frequency uj of an external a.c. electric field with polarization vector in the 
x-direction and wavelength large compared to £o- Right: The corresponding in phase 
charge response around the pinning site for frequency u) — 0.1 A. Gray corresponds 
to negative charge and black to positive charge. The arrows denote the total electric 
field vectors. The pinning center is a circular normal metallic inclusion with a radius 
£ p = 1.57£o, shown as the black circle. Impurity scattering is included in the Born limit 
with a mean free path of I = 10£o- The temperature is T — 0.3T C and the calculations 
are carried out in the high-K limit. 



In Fig. [T^ we show the total electric field (external plus induced) acting on 
the quasiparticles in units of the external field. For lu > 2 A the total field is 
approximately equal to the external field. However, at frequencies u> < 2A an 
out-of-phase response develops. For lu <C A, the total field in the pinning region 
approaches twice the external field, and the out-of-phase component vanishes. 
In the intermediate frequency region, O.IZi ^lu < Z\, both in-phase and out-of- 
phase components are comparable. The right panel of Fig. [l^ shows the charge 
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distribution for uj = O.lzA which oscillates out of phase with the external field. 
A dipolar charge distribution accumulates at the interface between the super- 
conductor and the normal inclusion, oscillating at the frequency of the external 
field. At the center of the pinning site the out-of-phase component of the field is 
nearly zero at low frequencies (see also the left panel) . The induced charge which 
accumulates is of order of eA/Ef within a region of order £q in each conducting 
layer. This charge is a factor of Ef/A larger than the static charge of a vortex 
that arises from particle- hole asymmetry 32 33J34|. 



4.2 Local Dynamical Conductivity 

Because the total electric field varies on the scale of a coherence length, the 
current response expressed in terms of the the total field is nonlocal in the 
intermediate clean regime. However, we can define a local conductivity tensor in 
terms of the response to the external field, provided the external field varies on 
a length scale large compared to the coherence length, 

J„{R,u) = <j^{R,u)E^{w) . (41) 




Fig. 13. Imaginary (left) and real part (right) of the local conductivity (multiplied 
with u> for convenience) in the center of an unpinned s-wave vortex (circles), and at dif- 
ferent distances from the center along the y-axis: l.l£o (squares), and 2.2£o (triangles), 
as a function of frequency u) of an external a.c. electric field with polarization vector 
in :r-direction and wavelength large compared to £o- Impurity scattering is taken into 
account in Born limit with a mean free path of I = 10£o- The temperature is T — 0.3T C . 
Calculations are done in the high-K limit. Also shown are the response for a homoge- 
neous s-wave superconductor (dotted, denoted 'BCS'), and the Drude conductivity of 
the normal metal (dashed). 

Figures [l3] and [jd] show results for the conductivity, mi, in the vortex core 
region as a function of frequency for both unpinned and pinned vortices. For the 
pinned vortex the radius of the pinning center is £ pin = 1.57£o- 
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First consider the unpinned vortex. The absorptive part of the conductivity 
(right panel of Fig. |l3|) is strongly enhanced at the vortex center compared to the 
normal-state Drude conductivity. The reactive response exhibits a maximum at 
a frequency determined by the impurity scattering rate. A few coherence lengths 
away from the vortex center the real part of the conductivity changes sign at 
low frequencies. This is a signature that energy is transported by vortex-core 
excitations away from the vortex center producing "hot spots" outside the core. 
The net dissipation is determined by inelastic scattering processes in the region 
around the vortex core. At distances of order a coherence length or so from the 
vortex center there is also structure in the conductivity spectrum at higher fre- 
quencies reflecting absorptive transitions between quasiparticle excitations with 
energies corresponding to the Van Hove singularities in the local density of states. 
The maxima in the absorptive part of the conductivity at y = 1.36£o and 2.72£ 
from the center correspond to the energy level separation between the Van Hove 
peaks above and below the Fermi level shown in Fig. [|. 




Fig. 14. Like Fig. ^ for a pinned vortex with a pinning radius £ p = 1.57£o- As 
pinning site a circular normal conducting inclusion located at the center of the vortex 
was assumed (see Fig. p^ ). 



The conductivity spectra for a vortex pinned by a metallic inclusion at its cen- 
ter is shown in Fig. [14[ The most significant difference compared to the unpinned 
vortex occurs at frequencies u> < 1/r. The absorptive part of the conductivity 
(right panel) is reduced compared to that of the unpinned case at low frequen- 
cies. The three broad peaks in the absorption spectrum correspond to scattering 
and dissipation within the zero-energy resonance (the dominant low-frequency 
peak), transitions between the zero-energy resonance and the continuum (the 
peak near to ~ 1.5A), and pair-breaking transitions from the negative energy 
to positive energy continuum states (broad peak above 2 A). The other notable 
feature is the reactive response at low frequencies which becomes negative in the 
low frequency limit, corresponding to superflow in the core that is counter to the 
induced supercurrent outside the vortex core and pinning center. This counter- 
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flow is required in order to satisfy the conductivity sum rule. The counterflow in 
the vortex center is also present for unpinned vortices, but at smaller frequen- 
cies for this particular impurity scattering rate. The low-frequency counterflow 
is also related to characteristic current patterns associated with low-frequency 
vortex dynamics which we discuss below. 




Fig. 15. a.c. current density pattern for an pinned s-wave vortex with a pinning radius 
£ p = 1.57£o- As pinning site a circular normal conducting inclusion located at the center 
of the vortex was assumed. The border of the pinning center (gray) is shown as black 
circle. The frequency of the external electric field with polarization vector in :r-direction 
is chosen lu — 0.1 A. Left picture for out of phase (reactive) response, right picture for 
in phase (absorptive) response. Impurity scattering is taken into account in Born limit 
with a mean free path of I = 10£o- The temperature is T = 0.3T C . Calculations are 
done in the high-K limit. 



4.3 Induced current density 

Results for the a.c. component of the current density near a pinned vortex are 
shown in Fig. [lB] for ui — 0.1 A In addition to the a.c. current there is the time- 
independent circulating supercurrent around the vortex center which adds to the 



current shown in Fig. 15. The current response shows a dipolar pattern, which is 
also observed for unpinned vortices. The in-phase current response (right panel) 
indicates a region of strong absorption within the pinning region (j||E cxt ), and 
emission (j ■ E cxt < 0) in the region roughly perpendicular to the direction of 
the applied field several coherence lengths away from the pinning center. Calcu- 
lations of the energy transport current show that energy absorbed in the core 
is transported away from the vortex center by the vortex core excitations in di- 
rections predominantly perpendicular to the applied field. The net absorption is 
ultimately determined by inelastic scattering and requires integrating the local 
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absorption and emission rate over the vortex array. Note that the long-range 
dipolar component does not contribute to the total dissipation. Far from the 
vortex core the current response is out of phase with the electric field and pre- 
dominantly a non-dissipative supercurrent. Also note that at low frequency we 
clearly observe the counterflowing supercurrent within the pinning center. 

4.4 Summary 

The electrodynamics of the vortex state in the intermediate-clean regime is non- 
local and largely determined by the response of the vortex-core states. Transi- 
tions involving the vortex-core states, and their coupling to the collective motion 
of the condensate requires dynamically self-consistent calculations of the order 
parameter, self energies, induced fields, excitation spectra and distribution func- 
tions. The results of these calculations provide new insight into the dynamics of 
vortex cores in conventional and unconventional superconductors. 
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